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ABSTRACT 


In order to characterize the propagation conditions 
aeoma KNOWN patns at VHF and S Band frequencies, transmission 
loss predictions are produced by computer methods. An attempt 
is made to define the standard atmospheric conditions along 
these paths through the presentation of the statistics for 


normal and super-refractive propagation conditions. 
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le NT RODMET TON 


Numerous methods of estimating tropospheric propagation 
path loss are currently found in the literature. These 
include both graphical and computer techniques, of which the 
graphical methods are most widely presented. In conducting 
this particular study a program dealing with tropospheric 
path loss prediction (TROPOPLOT) was chosen as a computer 
method of analysis. This program was originally published by 
ESSA in a technical report [Ref. 6] and was subsequently 
modified for use at the Naval Postgraduate School [Ref. 7]. 
Tests were conducted by Longley and Reasoner [Ref. 8], in 
which the computer-predicted results of the program were 
compared with empirical data taken over a number of different 
paths. The results showed that TROPOPLOT can provide 
reasonably accurate predictions of path loss within certain 
meniotraints. 

The paths selected for analysis in this study are located 
in the Puget Sound, Washington area. They comprise both 
existing and proposed communications links connecting the 
underwater range facilities operated by the Naval Torpedo 
station at Keyport, Washington. The range operation centers 
are currently located at Winchelsea Island, British Columbia 
and on Zeiatched Point near Dabob Bay. The present communica- 
tions requirements include on-range communications with both 
range control vessels and the submarines and aircraft conduc- 


mmig tests, the capability to monitor telemetry information, 


a2 





and the maintenance of telephone links with other sites. 
Proposals for expansion of the range facilities as the 

TRIDENT program progresses include possible at-sea range 

mem cies and digital data @eansmission to a central computer 
processing center located at Bangor or Keyport. These pro- 
posals would increase the bandwidth requirements of the links 
and thus necessitate a change from the present VHF frequency 
range to an S-band (2 GHz) line of sight link. 

The terminal and repeater sites are shown in Figure 1. 
iiemexisting VHF link consists of paths from Winchelsea to 
Lookout Mt. (repeater) to Bangor and from Zelatched Point to 
eager. the proposed link at 2 GHz runs from Makah to Striped 
Peak (repeater) to Mt. Constitution (repeater) to Gold Mt. 
(repeater) to Bangor or Keyport, as well as from Zelatched 
Point to Bangor or Keyport. Each of these paths 1s computer- 
analyzed in Section II. 

Maeaddition to computing path loss Jeune Wistorie¢d! 
meteorological data for this area was examined in an attempt 
Mmemenaracterize the "typical" propagation conditions found in 
this vicinity. This was accomplished by means of the 
Integrated Refractive Effects Prediction System (IREPS) which 
Provided output based on data taken from upper air sounding 
mratvons in the Puget Sound area. This information could also 


be useful in the determination of link service probability. 
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A listing of a computer program designed to plot path 
profiles and to aid in the determination of several input 
variables is included in the appendices, as is an explanation 
of the modifications required in order to use detailed terrain 


profile information in the TROPOPLOT program. 
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MiP CUE n=n DE DePALH LOSS PREDICT RON 


A. PROGRAM DESCRIPTION 

Pee UrL OT Ro nides as an output a measure of the long- 
term attenuation over a given path. The program was primarily 
designed for use in cases where detailed terrain information 
was not available but can be used accurately with user- 
supplied data after a slight modification. The required 
input parameters include frequency, antenna heights, antenna 
gains, line losses, receiver sensitivity, path length, 
Bgtdadrization, transmitter power out, surface refractivity, 
conductivity, permittivity, and a measure of the terrain 
roughness Ah. These parameters are for the most part wel] 
defined in Refs. 6 and 7 and are easily determined except in 
the case of user-supplied terrain profile information. For 
this sort of input there is some ambiguity concerning the 
selection of values for o, «, ah, and the effective antenna 
heights. 

In considering a path which includes both poor ground 
surfaces (c=4, g=.001) and sea water surface (c=81, o=5) the 
Bemticular value to use as input is not specified, although 
a completely oversea path requires an adjustment to the 
program not available in the version used for this study. 
Thus the values of e« and o chosen were based on the corres- 
ponding value for that portion of the terrain that ccnstituted 
tne largest portion of the dominant reflecting plane between 


ime transmitter and receiver. 
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The path geometry considered is shown in Figure 2 and 
indicates Nal ,2 as the structural transmitter (receiver) 
height above ground. This input parameter is converted into 
an effective height, 1h by the program. It should be 
noted that Longley and Rice consider two cases: one for 
random antenna siting in which structural heights, ee are 
considered equal to the effective heights, Net .23 and another 
in which the sites are carefully selected, as inthe case of 


radio relay links, and the effective heights are larger than 


the structural heights as shown below: 
2 + kexp(-2h 1 9/4h) meters 


The variable k is considered to have a maximum value of 
50 as determined by the author's study of varied terrain 
conditions. The method presented below for the computation 
of k is valid only for antenna heights less than or equal to 
10 meters. In Callaghan's version [Ref. 7] this determination 


WdaS: 


kK =] + Asin(ah , 4/10) FOr GOs hyl.2 on 
= otherwise 


Since most of the antennas in the links under study had 
menuctcural heights in the neighborhood of 50 meters, it jis 
eeerous that predictions based on strictly Nal ,2 would in this 
case, produce an erroneous (higher than normal) prediction of 
Path loss. This factor is recognized by the authors as an 


area requiring better definition and as a primary source of 
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Meedvetion error. For this reason the modifications to 
TROPOPLOT which are detailed in Appendix B were necessary. 
The effect of varied antenna heights in the input of the 
unmodified program 1s shown in the output data presented 
later in this section. 

The description of the terrain is accomplished by a 
statistical quantity, Ah, which is defined in Ref. 6 as "the 
asymptotic value of the interdecile range, Ah(d), of terrain 
heights above and below a straight line fitted to elevations 
above sea level. The parameter Ah(d) is calculated at fixed 
distances and its median value usually increases with path 
mengtn to Ah. For a single path profile these definitions 
are not adequate. The asymptotic value, Ah, used in this 
Study as determined by taking the interdecile range of the 
difference between a straight line fitted to the path profile 
points and the path profile points themselves. The results 
of this estimation of Ah correlate favorably with those 
values listed in Table 1 [Ref. 7] which shows estimated values 
of Ah for particular types of terrain. Figure 3 shows an 
example of the straight line fit to the terrain profile between 
Lookout Mt. and Bangor. The details of the computation are 
found in Appendix B. 

Having thus considered the variations in the input 
Parameters the program constraints should be noted. TROPOPLOT 


1s designed for use within the following constraints: 


eg) 
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Examples of Least Square Line Fit to Terrain 
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Parameter Range 


frequency (f) 20 to 40,000 MHz 
antenna height ee G2 5. 1to 3000 m 
astance (dist) 1 to 2000 km 
Sietace refractivity ce) 250" tO eeU0 Neunits 


ine antenna Siting 1s subject to the following conditions: 
ie: The angle of elevation, Bal .2? of each horizon 
ray from the horizontal should not exceed 12°. 
ee The distance from each antenna to its horizon 
(di yoo) Should not be less than 1/10 or more 
than 3 times the corresponding smooth earth 


distance (dic). 


These limitations can be ignored if the values of Ms a 
and B12 are entered directly into the program. As is the 
case with Ah and Na, 2 the method of accomplishing this is 
Someained in Appendix B. 

mavying resolved, or at least indentified, the potential 
ambiguities in the input parameters, it is useful to briefly 
describe the output format of the program. As shown in 
mise 4-7, the output of TROPOPLOT consists of two tables 
containing transmission loss versus distance and signal 
Semagth versus distance, a plot of the transmission loss 
information and a printout of the calculated and input values 
of many of the program parameters. Some of the quantities in 
this figure may require some explanation as to their meaning 


Or derivation. JThese include the variables AE, AES, AED, MS, 


wee nA, DX, Kl, KZ, and ALS. 
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Inthe process of describing these parameters a brief 
outline of how TROPOPLOT determines path loss is presented. 
A more detailed exposition can be found in Ref. 6. The 
median reference value of attenuation below free space, eee 


meomeompuced first. The reference value of transmission loss, 


L__, then becomes the sum of the free space attenuation, L 


EY bf? 


and the reference attenuation, ese 


Lor ” Lb F 7 Aor 08 


where the free space loss is defined as, 


L eo 45>, F20 logy, (f in MHz) + 20 log, ,(d in km). 


bf 


The reference attenuation, Pee is determined using one 
of three subroutines depending on the particular mechanism 
feemmopadation; LOS for line of sight modes, DIFF for 
iieniraction, and SCATT for tropo-scatter. Jwo ray optics is 
used to compute line of sight paths while diffraction paths 
are assumed to be over a double knife edge and a value of 
diffraction attenuation below free space, AG? is computed. 
In the case of scattering the attenuation variable is designa- 
ted as A: The reference attenuation ees is determined by 
the smaller value of AG and A. In the 1 ine~or svahteease 
the attenuation is calculated for 2 values of distance 
(do sd.) for which line of sight propagation is valid to 


produce corresponding values of attenuation aN and Ay. The 


diffraction attenuation, ine 1s computed at the distance 
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dy. and together with Ao and A, is used to determine the 


eehOOe S , Ky and Ko: of a smooth curve of the reference attenua- 


eine versus distance over the range | < d < od 


pie = a i k, (d-d,) co Ko log, ,(d/d.) aise, 


The output parameter Ean 1s defined as 


A =A_- k.d_ - Ks 10g, 94 (i.e., d=1) 


e 0 are o 
d 
moepor | < d < a 
er = A és kd i Ko Tog, 94 


When the diffraction attenuation is computed the output 


variables ae and m, come into consideration. The diffraction 


d 
attenuation is computed as the weighted average of the 


diffraction attenuation over smooth earth, A? and the attenu- 


ation over a double knife edge surface, Avs where 


A = (1-w)A, a WA. dB 


A description of the method of weighting these estimators~ 


aaceecal culating A and Ay 1s contained in Annex 3 to Ref. 6. 


iemaitfraction attenuation, A 1s determined at two 


d 3 
A? im tne: tar diffraction regqiom, sand a 
straight line through the points (A 


@rstances, ue aad d 


3°43) and (A, sd,) 1s 


defined as 


ma = (A, - Aa)/(d, - d,) dB/km (slope) 
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and 


Ba S Ac, + Ay - mad, (intercept) 


where A, is a “clutter factor" (< 15 dB). 


If the scatter attenuation is less than the diffraction 
attenuation, as is sometimes the case in trans-horizon paths 
where the distanced or the angular distance 9 is large, then 


2. i = Bie When the product of d in kilometers and 9 in 


maa@ians iS greater than 0.5, A. 1s computed at large distances, 


d. and de, and a straight line through the points (A d. ) and 


(Ag sd) ismderined as follows: 


5 > 


and 


A =A_=A + m_d d F2sha 
S ae 


The quantity oe is defined as the distance where the scatter 
Peeeomuation 1s equal to the diffraction attenuation. 

It can be seen from the above description and Ref. 6 that 
the actual determination of the output quantities is somewhat 
complex and for most applications the information contained 
in the graphs and tables describing path loss and signa] 
Strength is sufficient. Consequently the data presented in 
the next sub-section will primarily consist of derivations 


maom those portions of the output. Jf the quantities in 
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meeo «are important to tne user it should be noted that for 
certain input combinations the format statements contained 
in the TROPOPLOT program will produce an all asterisk, 

meee =X". printout indicating that the format statement 


governing that particular variable requires modification. 


B ACTUAL PATH LOSS MEASUREMENTS 
le Procedure 

In order to obtain the proper input parameters for 
entry into TROPOPLOT certain preliminary measurements and 
calculations were required. The first step was to obtain 
Saarts Of the path area! and from these determine the terrain 
profile over a great circle path between transmitter and 
receiver. The great circle path was approximated by a rhumb 
line for path lengths less than 70 kilometers and the avail- 
ability of 7%' topographic charts for the longer paths enabled 
the use of straight lines over the chart area. This profile 
MarecmatiOn was plotted on a curved earth's surface using 
the plotting program described in Appendix A. The value of 
the surface refractivity used in determining the effective 


earth's radius was the same as that used in the input to the 


IT Obtain aeecatdlogue for U.S. Charts write to 


Venver Diseri bution, seetion 

U.S. Geological Survey 

Denver Federal Center, Building 41] 
Denver, Colorado 80255 
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meogram and can readily be obtained from Ref. 9 or from 
mecual meteorological data. [f the charts of Bean et al. 


(1960) are used this value must be converted from N, to N. by 


the equation 


ae i exp(-0.1057 i) 


The value of h. used depends on several factors most important 
OF which is the mode of propagation. For most of the paths 
considered herein the value of nh. chosen was determined by the 
elevation of the lowest antenna. For trans-horizon paths a 
mean N. 1S computed using heights at the obstacle horizons. 
As was previously noted the values of Ne FOr 1nDUT Can wane 
meonmecoU to 400 N-units, so it 1s obvious that predictions of 
path loss for anomalous conditions are not readily obtained 
emioeuROrOPLOT which is designed for long-term and median 
input parameters. The effects of non-standard atmospheric 
conditions along these paths are considered in Section III. 

Apter plotting the path profile the required inpud 
Darameters were determined using the methods described in 
Appendix B and a number of parameters were varied to produce 
the results described below. 

a. Res Ul tse : 

The output of the path profile/routine for each 
path considered is shown in Figs. 8-19 while the corresponding 
input parameters used are contained in Table I. All the links 


were analyzed for path loss at both VHF and S-band to 


SL 





HEIGHT ABOVE SEA LEVEL ( m ) 


140 


120 


100 


OO 
oO 


iZelatched Point (60.96) 


o 
O 


p 
© 


20 


3 Keyport: | eed 1) 


5 10 15 
DISTANCE (km) 


Figure 8 


Terrain Profile for the Zelatched Point to Keyport Path 


by 





HEIGHT ABOVE SEA LEVEL (m) 


Makah (365.76) 


300 


200 


100 


O 


Striped Peak 
(255.4) 


20 40 60 
DISTANCE (km) 


Figure 9 
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Terrain Profile for the Bangor to Gold Mt. Path 
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Terrain Profile for the Lookout Mt. to NOEF, 
Bangor Path 
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illustrate the increase in path loss when frequency is 
increased. Tne results are shown in Figs. 20-31. The quanti- 
ties (referred to in an above sub-section) which required 
interpretation as to the precise value to be used in detailed 
beiemainm Studies are varied for a single path so that the 
effects of this variation can be shown. A sample path 
(Lookout Mt. to Keyport) was selected for this purpose and 
the effects of several different antenna heights in an 
unmodified version of TROPOPLOT are shown in Fig. 32. These 
Same antenna heights are then entered into the modified 
eoodmam and the results shown in Fig. 33. To show the effects 
of variations in conductivity and permittivity o and ¢« are 
varied over the same path with the results shown in Fig. 34. 
The effect of variations in the value of Ah is shown in 
migmmeoo. ihe quantity Ah was computed by taking the mean 
value of Ah over all the paths considered, An, is the tabulated 
Sramicic10On Of terrain irregularity found in both Ref. 6 and 
mer. 7, An. is the parameter computed by the methods of 
Appendix B. This comparison is intended to show the degree of 
precision required in choosing a parameter when uncertainty 
Over tne means of selection is present. 

The importance of detailed path profile information 
is pointed up in the Zelatched Pt. to Keyport path (Fig. 8). 
If TROPOPLOT is used without modification, errors in prediction 
can occur at 2000 MHz the results of which are shown in 


megs 20. Note that this output predicts free space attenuation, 
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Figure 25 


Transmission Loss vs. Distance for the Striped Peak to 
Mi. Constr tievon Path 
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Transmission Loss vs. Distance for the Bangor to 


Gold Mt. Path 
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Transmission Loss vs. Distance for the Keyport 
to Gold Mt. Path 
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Figure 28 


Transmission Loss vs. Distance for the Lookout Mt. 
to NOEF, Bangor Path 
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Figure 29 


Transmission Loss vs. Distance for the NOEF, 
Bangor to Gold Mt. Path 
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Figure 30 


Transmission Loss vs. Distance for the Zelatched Point 
to NOEF, Bangor Path 
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Figure 31 


Transmission Loss vs. Distance the Lookout ME. 
to Bangor Path 


56 








Debs ie Oem. 
ek ) 


180 ie 


160 a H1,2=ht. above 
“ Sea level 
o ww 
; re 
ay re C) 
a a _ Tabu) & Values 
= 140 / y Co 
Y Co & 
2 aC eS 
zs Pc es 
= J 
= 120 ao 
F eS 
se jff 
D yy 
Ee Ow 
= 100 
80 
25 50 75 LOO as 


DISTANCE (km) 
Figure 32 


Antenna Height Variation for an Unmodified Version of TROPOPLOT 
Transmission Loss vs. Distance 
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Antenna Height Variation for the Modified Version of TROPOPLOT 
Transmission Loss vs. Distance 
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Variation in € and o for the Lookout to Keyport Path 
Transmission Loss vs. Distance 
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Variation in Ah for the Lookout Mt. to Keyport Path 
Transmission Loss vs. Distance 
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when there is clearly a large attenuation of the signal by 
Memenoostacies. The same error can occur when the modified 
program is used; primarily due to the fact that the terrain 
is represented statistically. Large deviations from the 
range of terrain Neights occurring in mid-path may not be 
considered as an obstacle by the program when in fact it is 
pemeyor factor in the attenuation process. These smoothing 
effects can be alleviated somewhat by examining the terrain 
profile and performing graphical analysis, such as is 
mescrirped in Refs. 1-5, to check the effectiveness of the 
computer routine when its output is questionable. Another 
possible source of error can occur in short paths where only 
a single knife edge obstacle is present. In this case the 
predicted value of path loss is too high since a non-LQS path 
loss is calculated for a double knife edge case. 
Se Cone ius 1Ons 

HE is obvious that the quality of the outputvot 
this tropospheric propagation prediction program depends on 
which of several possible assumptions were made concerning 
the choice of input parameters. Clearly, if information is 
desired concerning a particular path or path area a path 
profile should be drawn and the ray path checked to insure 
that the values of the corresponding angles of elevation 
(depression) fall within the program constraints. As 
Omeviously stated, if these values exceed the limits then 
TROPOPLOT must be modified and a detailed terrain profile 


obtained. 
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ime crrect Of Tie Variations Of several of the 
parameters shown above allow some inferences to be drawn 
concerning the care with which the input variables must be 
chosen. It would seem from the path studied that variations 
Gao ana ¢ have the least effect on path loss (<1 db), varia- 
meomowim An nave only a small effect, particularly for longer 
distances, and variations in antenna height have the greatest 
effect on path loss. The unmodified version of the program 
appears more susceptible to these changes than the version 
modified to accept detailed profile information. The sensi- 
tivity of the transmission loss to antenna height value 
Variation was alluded to in an earlier section and thus the 
Suggestion that care be exercised in selecting these values 
is well taken if any degree of precision is desired. The 
detailed computation ofdah, on the other hand, is not necessary. 
By using the tabulated values corresponding to a particular 
type of terrain [Ref. 6] similar results can be obtained. 

The transmission loss predictions seem to be 
Mectrate, or at least within an expected range of values except 
monesmort diffraction paths. In order to assess the goodness 
Sieecwe orediction, comparison with empirical data for these 


paths should be accomplished in a later study. : 
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III. PROPAGATION IN A NON-STANDARD ATMOSPHERE 


The computer analysis conducted in the previous section 
produced a long-term median attenuation value as output and 
consequently required that standard atmospheric conditions 
be assumed. In providing input under this assumption, the 
"standard" atmosphere for the locations under consideration 
was characterized by the median value of minimum monthly 
wimtace retractivity at sea level (No) corrected to emitter 
elevation. This element of the prediction routine could, in 
some cases, lead to an over-optimistic expectation of link 
performance, particularly in a locale where anomalous condi- 
tions frequently occur. One method of avoiding this potential 
moueece Or error is to examine the statistical occurrence of 
these conditions and the severities of the effect of the 
anomalies on the link in question. For the purposes of this 
study the non-standard conditions dealt with involved the 
Super-refractive, subrefractive and ducting cases. The 
statistics used in this determination were derived from IREPS 
Mmeput and the effects of non-standard refractivity were 
modeled by assuming certain values of effective earth's radius 
for corresponding conditions. Prior to the presentation of 
these results the theory involved in the anomalous propagation 


problem is reviewed. 
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A. THEORY 

The primary medium through which electro-magnetic waves 
with wavelengths less than 1 or 2 meters propagate is the 
troposphere. This region, bounded below by the earth's 
Surface and above by the tropopause, 1s approximately 10 km 
thick and 1S characterized by a general decrease in tempera- 
ture with height up to a zone of constant temperature called 
miemtropopause. The tropopause 18 not a static boundary, but 
has a height which is variable with both time and latitude. 

The troposphere is usually assumed to be a lossless 
mpetiectric with »p = 1 and o = 0, The index of refraction, n, 


is then 
“ens 
At the earth's surface n has been found to equal approximately 


meoeo3. Since air with a higher water vapor content has a 


larger value of Herth! Helv ley wud e.) Ws. 


fewet air — °dry air 
and since the lower atmosphere usually has the higher water 
vapor content, the permittivity of the atmosphere 2x hip esac 
decrease with height to a value of unity. This effect causes 
the refractive index to similarly decrease with height. It 
has been ascertained [Ref. 3] that for temperate climates the 
average variation near the ground is 


i = 0.039x107° per meter 
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In order to avoid the use of such small numerical quantities 
meeceiable, N, designated as the refractivity or co-index of 
refraction, has been defined as 


6 


N = (n-1)x10 N-units 


wicesO within one kilometer of the earth's surface 


os = -39 N-units/meter 


Meotanddard atmosphere is then defined in Refs. 11 and 12 as 


mo | 1OWS: 
] eels 
n= | 
2 a 
eran Is m 
77 6 as p = pressure (mb) 
aa (p+4810 +) T = absolute temp. (°K) 


vapor pressure (mb) 


ByeassuUming a linear refractivity gradient it is possible 
to define an effective earth radius as was done in an earlier 
Sarc Of this report. This allows the rays representing the 
radio wave to be drawn as straight lines over a curved earth. 
An alternative representation which enables the effect of 
several different values of effective earth radius to be 
Shown on the same profile is to snow the rays as curves from 
transmitter to receiver over a plane earth. Details of 
constructing this sort of representation are presented in 


mets. 3 and 13. 
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Non-standard propagation can occur when the above condi- 
tions are not met. This is manifested by the manner in which 
the radio wave path is curved. nae momma! C@iNcd du) Ons seme 
decrease of the refractive index with height causes a downward 
bending of the wave path. Changes in the refractive index 
Super-refraction or in certain cases an upward bending of the 
ray which is termed subrefraction. If the refractive effects 
are severe enough, the formation of ducts occurs in which the 
majority of radio wave energy is trapped within a narrow 
megion. 

Ducting can occur as a result of a number of different 
atmospheric conditions. The primary requirement in producing 
this phenomena is that the curvature of the radio wave must be 
greater than that of the earth. This is produced by a rapid 
change in N with height which is caused by a sharp decrease in 
moisture with height (abnormal moisture lapse rate) and/or a 
sharp increase in temperature with height (temperature 
inversion). The formation of a ground- or surface-based duct 
can occur when warm dry air flows from land over water. 
Moisture evaporates from the water into the lower layers of 
meee oo!ling the air and producing an increased moisture lapse 
rate and a temperature inversion. The downward motion of warm 
dry air is frequently associated with the clear weather found 
on the eastern side of vast high pressure regions in the lower 
and middle atmosphere. Subsidence such as this can also 


meswlt in elevated ducts which are found to be strongest and 
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lowest in fair weather and highest and weakest near storms. 
Surface ducts are also formed in certain regions by nocturnal 
cooling but unlike the evaporation duct this sort of anomaly 
ienound over land, particularly desert regions. 

The opposite effect to super-refractive ducting is the 
upward bending of radio waves caused by an increase in 
BeteeaetiVity with height. This condition typically occurs 
foeemeist and cloudy regions. 

litemcetermining factor in all the above cases is the 
refractivity gradient. Table II below, found in Ref. 11, shows 
the range of refractivity gradients over which anomalous 
Pomagitions orevail. 


Ia te ii, Range of refractive index gradient for differing 
types of propagation conditions. 


Sub - Super- 
memo) tC Km) | Refractive | | Unstratified | Refractive 


| 1000-850 0-1.46 | -dn/dh<0 | 20<-dn/dh<60 | 100<-dn/dh 
‘| 850-700 | 1.46-3.01 | -dn/dh<0 | 20<-dn/dh<50 80<-dn/dh 
700-600 | 3.01-4.20 | -dh/dh<O | 20<-dn/dh<40 70<-dn/dh 


600-500 | 4.20-5.57 | -dh/dh<O | 20<-dn/dh<30 50<-dn/dh 


| 500-400 Dry =f lie -dn/dh<0 2Uec= di) dinae 40<-dn/dh 


| 
| The refractivity gradient can be related to the 
















effective earth radius and the resulting ray trace when 
| drawn on the terrain profile corresponding to this radius, 


will show the effect of the anomaly on a particular path. 
| 
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imemiol|lowing holds for most conditions [Ref. 3]: 


] 


erreceive fTadius =k pO ome aaa REGS a 
leeO Oe a Neaih 


meomecnis relationship it can be seen that for the standard 
refractivity gradient, which was defined as -39 N-units/km, 
the effective earth radius would be equal to 4/3. The 
presence of anomalous conditions thus cause an effective 
eartn radius greater or less than 4/3. The interrelationship 
of the refractivity gradient, the radius of curvature of the 
wavefront, the effective radius, and the assoicated types of 
Mearaction are shown in Table III [Ref. 3]. 

Note that there is some difference in the exact gradient 
value which defines super-refraction. The super-refractivity 
Meaditent referred to in Table II includes both the extended 
Gomge and the ducting conditions, while that referred to in 
mippe Ill includes only the ducting case. 

In the actual paths under study several values of 
meaective earth radius were used to simulate the effect of 
anomalous conditions. Table IV [Ref. 15] provides a guide 
to the use of k as an estimate of propagation conditions for 


99.9%-99.99% path reliability. 
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Lookout Mt., and Bangor to Lookout Mt. are shown for two 
memes Of K: 0.5 which represents the difficult case as either 
when one antenna is inside a duct and the other outside, or 
ier steretractive condition, and -3.5/7 which represents a 
subrefractive gradient of -200 N-units/km. The effect of the 
change in earth radius shown in Figs. 36-41 is most noticeable 
in the Keyport to Lookout Mt. path. Using graphical methods 
Seanelysis found in Refs. 1 and 3 the approximate degradation 
Meilink performance can be estimated. While little change is 
Mimserved in the Keyport to Gold Mt. link due to the snortness 
of the path length, the Keyport to Lookout Mt. path with K = % 


Bea bits an additional loss of up to 60 db over that found in 
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Figure 36 
| Path Profile for Lookout Mt. to Keyport (K=1/2) 
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Figure 37 


Path Profile for Lookout Mt. to Keyport (K=1) 
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| Path Profile for Keyport to Gold Mt. (=) 
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Figure 39 
Path Profile for Keyport to Gold Mt. (K=1) 
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| 





Sandard atmospheric conditions. This same path however 

suffers no apparent degradation for K = 1 as the Fresnel zone 

clearance remains greater than .6. The Lookout Mt. to 

Bangor path degradation was found to be 12-30 db for K = & 

mimo-ho db for K = |. From these results it can be seen 

that changes in effective earth radius can have a significant 

effect on path loss, particularly as the path length increases 
The profile representation using N data is the most 

commonly used in the Jiterature, however the depiction of 

meena or trapping types of profiles is often done with 


modified or M-unit profiles where 


actual earth radius 


a 
m= ny + 442 x 108 


= 
tt 


height above sea level 


This adjustment to the N-profile is used in the anomalous case 
Since straight rays above a curved earth now become curved 
rays above a planar earth and the duct can be more clearly 
wiowmee ine variability of M curves corresponding to a 
particular propagation condition is shown in Table V which 

was derived from Ref. 10. 

The M-curve is thus a transformation in which tne 
relative curvature between the normal of the electro-magnetic 
wavefront and the surface of the earth is unchanged. A duct 
is formed whenever the M curve has a relative minimum. The 


base of an elevated duct is the height on the M curve below 


the minimum that has the same M value as the relative minimum. 


Vi 





meremo SUCNh height exists, the base of the duct is at the 
Geodiaeand a so0-called ground based duct exists [Ref. 12]. 
The main disadvantage in the use of M data is that above the 
duct N is grossly overcorrected. 


Table V. Variability of M curves with type of transmission 
eOmadi tions . 


Type of 
M M-Curve [as Giessen 
Increases linearly SbCalgne linvewnw ion Stamaard 
positive slope 
Slope decreases near Substandard 
eartn surface with 
rays -euny ingeup 
Slope increases near Superstandard 
Surface with less 
upward curve 
Vertical, no curve Very great 
coverage 
Decreases with Slope linear and Trapping or 
heignt negative rays duet 1G 


curve down 





“ Having demonstrated the effect of anomalous conditions 
on several propagation paths and considered the theory 
involved the next step is to determine tne frequency with 
which these conditions occur. Tne method chosen to accomplish 
this was to access the [REPS historical file for the geograph- 


ical area in question as described in the following section. 
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a Ime GRATED REFRACTIVE EFFECTS PREDICTION SYSTEM fiREpS 
[REPS iS an experimental system ultimately designed to 
provide the capability for on-board assessment of the effect 
of atmospheric anomalies on sensor performance. The hardware 
portion of the system consists of a mini-computer with an 
interactive graphics terminal on which the required inputs are 
entered and the results displayed in an easily understood 
graphic format. The information stored in the systems memory 
includes a refractivity library containing long-term meteoro- 
Hogdiecdie statistics on the occurrence of ducting conditions as 
[eriernon of location. This library data is obtained from 
radiosonde stations throughout the world and can be augmented 
by on-scene refractometer or radiosonde data input in order 
momedadllor the output to a particular location. In the case 
of the Washington paths, there were only two upper air sounding 
measepons in the vicinity of the sites. Qne is located at 
Tatoosh Island (48°-24'N, 124°-42'W) and the other at 


Quililayute (48°-00'N, 124°-36'W). The historical propagation 


/Msummaries for each location are shown in Figs. 42 and 43. 


Note that while the statistics generated are geared toward 
madar performance, the frequencies to be used in the link are 
included so that some information concerning potential 
emomalies can be derived. The historical meteorological data 
momeeach station, shown in Figs. 44 and 45, provides informa- 


micmmeconcerning median conditions in the vicinity of tne 


er or qedetai led description of thevsystem see Rep. 10 


Ao 





—————————— 


sounding station and shows the seasonal and diurnal variation 
foeemis data. Based on this data, the program then can 
Memerate a series Of profiles. Figures 46 and 4/7 are M and 
feororiies in the vicinity of the two sounding stations. 

Also available in the program output capabilities is a ray 
Mmeaece diagram from a particular emitter. This is useful in 

a qualitative sense in that it shows potential areas of non- 
coverage which may be important in maintaining a ground to air 
communications link. The first set of traces demonstrates the 
menect on the ray trace diagram of selected sites for the 
median ducts found in the Tatoosh Island profile and are 

mown in Figs. 48-52. A second set of traces are shown in 
Figs. 54-57 exhibiting the effect of an elevated duct between 
mieO0 and 1400 feet. The profile responsible for this condi- 
mon 1s shown in Fig. 53. The third set of traces show a 

duct between 2300 and 2900 feet in Figs. 59 and 60 and the 
merresponding profile in Fig. 58. 

Besides the above mentioned output format the IREPS 
program can also provide path loss information, with a 
Metection threshold for a particular sensor superimposed on 
maeedistance vs. dB display and also a coverage diagram for 
a particular sensor. These two output diagrams were not 
produced for this study but an example of pach wseincluGded sin 


Figs. 61 and 62 [Ref. 10]. 
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Mmis system is not yet Operational anid isisstq isla en— 
going testing. The value of the output formats in evaluating 
link anomalous propagation potential is clearly obvious. If 
data were available for each of the sites analyzed in the 
study a more exact picture of the statistics of propagation 
would be available. 

Based on the available data for the Washington area 
several inferences can be made concerning anomalous propaga- 
tion. For the frequencies in question the probability of 
extended ranges or ducting is between 8 and 10 percent while 
no statistics are presented for below normal refraction. For 
the locations from which the soundings are available it can 
be seen that surface ducting occurs most often in the daytime 
during Spring and that the ducts are predominantly less than 
12 meters in thickness. The lowest occurrence of surface 
ducts at Quillayute is in the spring night which also has the 
highest incidence of elevated ducts. At Tatoosh which is a 
more classically marine climate, there are no ducts observed 
to occur during the autumn and winter either at night or 
during the day. The differences in the median refractivity 
meesea level are not significant and both values, 332 for 
Tatoosh and 324 for Quillayute, would produce an effective 
earth radius greater than 4/3, and consequently better than 
average propagation conditions prevail. Until meteorological 


data is available for the actual sites the ray trace diagrams 


Oy 








can provide only an indication of the effect of a certain 
profile on the emitter in question rather than actual informa- 


tion concerning real conditions present in the link components. 
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APPENDIX A 


Path Plotting Program 

Iigeorder to plot a path profile on a curved Garth surface 
either commercially available profile paper must be used or 
the profile can be manually or computer generated. Instruc- 
tions for the construction of path profile paper for various 
values of earth radius can be found in Ref. 13. Computer 
generated profiles can be easily produced if a plotting 
package is available. For the purposes of this study a 
subroutine known as DRAWP was locally available. Qther 
alternative plotting methods included the TEKTRONIX 4012 and 
HP 9830 plotting packages, both of which are available at the 
Naval Postgraduate School. 

The program used made use of a transformation which 
allowed the terrain profile to be plotted on a linear graph 
by adjusting the height at a particular distance for the 
effective curvature of the earth. The profile is drawn by 
the CALCOMP plotter along with the earth surface at sea level. 
The program also fits a line to the terrain data by the method 
of least squares and forms the difference between the value 
of height found on the line and the actual value of the 
terrain height at the same distance. Aside from the actual 
terrain data points the only input parameters necessary EO 
obtain this output are either the refractivity gradient or 


the surface refractivity, and the number of points input. 
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The output provided by the program includes a plot of the 
terrain on the effective earth surface, a table of height vs. 
ieee distance, a 11st of the terrain heights in kilometers, 
the effective earth radius, and the difference between each 
point on a least square curve and the actual terrain data 
points without the correction for earth curvature applied. 
In order to perform the curve to fit terrain referred to in 
Appendix B, a new set of points must be used taking only 
those the points visible to both transmitter and receiver instead of 
the entire terrain profile. 

The format of the input data deck is presented SeTONE 
piece Data Card: 


Column 1: Code "1" surface refractivity will be 
used as input 


Code "0" refractivity gradient will be 
used as input | 
Second Data Card: 


Column 1-10: Number of data points 
(right justified integer < 200) 


Column 10-20: Refractivity or gradient 
(right justified integer) 


Third Through Last Data Card: 


The distance points in kilometers are listed in order 
using real numbers (i.e., F10.5 format) in columns 
Hee Wil 20 ema nee Ole7 Ow Ae heights above sea level 
«1 feet follow directly after the distance listing 

an the same format. A new card is not started unless 
the last distance point was listed in Columns 61-70 


An example of the latter cards is shown below 


LO5 





For the table of distance vs. height below 


Distance Height 
(km) (feet) 
On0 420. 
0) 500. 
2.0 oe 
6.0 125.5 
4.0 780. 
5 A10, 1000. 


The final data cards would be as follows: 


Column: 1-10 11-20 21-30 31-40 41-50 51-60 61-70 


Cong. 3 0.0 ie 2.0 a0 4.0 5d 420. 
Card 4 500. 10: V2Z5e5 2/7 60.> = 1000: 
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APPENDIX B 


INPUT PARAMETER COMPUTATION 


A. DETERMINATION OF DELTA-H 

The program described in Appendix A was used to produce 
from the path geometry the parameters diy 2: 8 o1,2 and Ah. 
The computation of Ah is only partially computer-aided by the 
miming Of a straight line to the uncorrected profile heights 
and subtracting the height of the line from the corresponding 
berrain height. This difference is output in the array D{1). 
The interdecile range is then extracted manually from this 
data by finding the value of the lower boundary of the 90th 
percentile and subtracting from it the value of the upper 
boundary of the 10th percentile. An example of the fit of 
this curve to actual terrain data is shown in the main body of 


mimes study. 


B. DETERMINATION OF HE1,2 
For line of sight paths the following determination was 
made for effective antenna heights: 


ie A straight line is fitted to the corrected 
terrain profile heights. 


Oe. The line is defined by 
B(I) = AO + Al(X(1I)) 


in the program defined in Appendix A and should 
be modified for earth curvature by 


B'(1) = B(I) - X°/2A 





oe The points chosen in the actual profile to which 
this line is fitted should be only those points 
which are visible to both transmitter and receiver. 

4. The effective antenna height is then the difference 

between the fitted line and the height of the 
antenna above sea level. 
The procedure described above defines a reflecting plane 
between the transmitter and receiver. 

If, on inspection of the actual terrain and the least 
Square curve, a good fit does not appear to have been obtained 
due to the roughness of the terrain other methods of obtaining 
effective antenna heights are available. If the foreground 
of the antenna represents a good reflecting surface then the 
antenna height over ground can be used, or multiple least 
Square curves can be used if the path consists of a number of 
reflecting planes. In the event that terrain characteristics 
defy this sort of analysis, curves and formulae are available 
in Annex III to Ref. 5 which provide alternate methods for 
computing effective height. 

The curve to fit terrain is also used in the determination 
of effective heights for the case of knife edge diffraction 
which is assumed to occur when the path length is less than 
the smooth-earth radio horizon for each antenna. In this case 
two curves are fitted to each of the diffracted paths forming 
a reflecting plane for each antenna. 

For trans-horizon paths the effective heights are obtained 


by considering the actual height of the antenna above sea level 


anid subtracting from that height the average height of the 








terrain between the transmitter and receiver. To lll us trate 
this, consider a link with transmitter height above sea level 
— The heights at N equi-distant points are selected and 
the mean of the central 80% of these values is computed to 


produce hy where 


(Oo 
<— 


ah ee sa 
io > een ti 
l= 


i 
cal 
as 


27 
Evian = £8 22 
eee = 31; fp = oe DD hy, (2 0p eee ,30 
i=3 


The effective transmitter height is then 


h = h - fh Or Beech 


e | to t it to 


If the mean value is greater than the height above sea 
level then the structural height of the antenna above ground 


1s used. 


CZ DETERMINATION OF DL1,2 and TE1,2 

ime angular distance @ is readily obtained from) the 
geometry of the profile as shown in Fig. 2. 961? and 862 
are measured and 7 is calculated in the program as the 
maximum of either Oa] 7 849 oF -d, /a. The distance dy is the 


sum of the distances to the obstacle horizons, dy and di 5 


209 


Di MODIFICATION TO TROPOPLOT 


With the above means of determining the parameters 


eae” 


i: 


d 


10 
a 


ES 
14 
wo 


lye 


TROPOPLOT was modified as follows: 


Delete the following statements from the main 
pmo died ii: 

eer hoe BE eo) ) > GOs dia 

Me eCHiNG.GE «2 s) 2ORe HIG yEEsbe)) eGOn io ale 

Lips 97.0 

GON Toe 10 

Ale = OF DS IN (3. 1415927461 Gale.) 

H1E = H1G+Z1*DEXP(-2.0*H1G/DH) 

PeeHi2G ale ec sO) Woon 10s Io 

be Hee. Ge.2.0.0R H2G.LE.5-0)) GOT Omi 

ee: 

GO TO 14 

Z2 = 1.0+DSIN(3.1415927*H2G/10.0) 

H2E = H2G+Z2*DEXP(-2.0*H2G/DH) 

CONTINUE 

DLI DLST*DEXP(-.07*DSQRT(DH/DMAX1(5.0,HI1E) )) 
DL2 DLS2*DEXP(-.07*DSQRT(DH/DMAX1(5.0,H2E) )) 
TEI (.00065/DLS1)*((DLS1/DL1-1.)*DH-3.077*HIE) 
nec (.00065/DLS2)*((DLS2/DL2-1.)*DH-3.077*H2E) 


Insert the following statements 


909 FORMAT (4F10.5) 

PERO ee SOS UDI Me sete tren tee 
after the statement 
22 FORMAT (4F10.5). 


Having accomplished these insertions and deletions the input 


card deck 


as shown in Ref. 7 is changed as shown below for 


user-supplied data only: 





Popes t Data Card: 


Column 


01-10 


Code: 4 


Column 


11-20 


Code: 0O 


Column 


Second Data 


Column 


Column 


Column 


Column 


Column 


Column 


21-30 


Garda: 


01-10 


11-20 


21-30 


31-40 


41-50 


51-60 


“Type of Terrain" - Integer 


User Terrain Data Used 
No data suppression 


User Terrain Data Used 
No "output parameters" printed 


"DB Loss" - Integer 


No DB data desired 
(omit 2nd data card) 


DB data desired 


"Distance Between Antennas" - Real 
(1-2000) KM 


"TX Power Out" (Watts) 
Positive Real Number 


"TX Antenna Gain" (db) 
Positive Real Number 


"RX Antenna Gain" (db) 
Positive Real Number 


"Transmitter Line Loss" (db) 
Positive Real Number 


"Required Receiver Line Loss" (db) 
Positive Real Number 


"Raceiver Sensitivity" (dbm) 
Positive Real Number 








ied Data Card: 
Column 01-10 


Column 11-20 


Column 21-30 


Column 31-40 


Fourth Data Card: 


Column 01-10 


Code: 01.00 


Code: -1.00 


Column 11-20 


Column 11-20 


Column 31-40 


Eievenmepata Card: 
Column OI-10 


Column 11-20 


Column 21-30 


Column 31-40 


/SUnTdAce Re TYractivi tye 
Integer (250-400) N-units 


“Surface Conductivity" 
Real (mhos/meter) 


"Relative Dielectric Constant" 
Real 


"Interdecile Range" 
Delta-H 


"Antenna Polarization" 
Real 


Vertical Polarization 
HORizontalvurpelanizauion 


"Frequency" Real 
(20-40 ,000) 


"Calculated Antenna Height of 
Transmitting Antenna" 
Real (0.5-3000.0) Meters 


"Calculated Antenna Height 
of Receiving Antenna" 
Real (0.5-3000.0) Meters 


"Transmitter Obstacle Horizon, DL1" (km) 
Real 


"Receiver Obstacle Horizon, DL2" (km) 
Rea | 


"Transmitter Elevation (Depression) 
Angle" (radians) 
Real Positive (Negative) 


"Receiver Elevation (Depression) 
Angle TE2" (radians) 
Real Positive (Negative) 
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